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ABSTRACT.- Rapid and accurate identifications are crucial for biodiversity assessment. Yet, traditional methods 
for species identification have some limitations. In this study, we tested the efficiency of mitochondrial 
cytochrome c oxidase I barcoding sequences for species identification and documenting diversity of aquatic 
insects in the rice fields of Thailand. Considerable success rate (8096) for species identification was found among 
the species of the order Odonata. Unidentifiable specimens of immatures were successfully associated with 
conspecific adults or by matching with reference sequences in the public DNA barcoding library. However, some 
specimens were ambiguous, possibly due to incomplete lineage sorting of closely related species or erroneous 
identification of the sequences in the public database. The technique was less successful for other insect orders 
because a lack of reference sequences in the DNA barcode library limits the utility of DNA barcoding. The 
Poisson tree process and Automatic Barcode Gap Discovery species delimitations revealed that the number of 
species recognized is more than twice that based on morphological identification. Therefore, DNA barcoding has 
potential for use in species identification and biodiversity assessment of the aquatic insects in the rice field 
ecosystem. 
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other inappropriate managements could 
reduce biodiversity in agricultural ecosystems 
(Wilson et al., 2008; Watanabe et al., 2013). 
Thus, information regarding species diversity 
in the rice field is needed for monitoring as 
well as appropriate management. 

Aquatic insects are among the major 
components of the rice field faunas and play 
diverse roles in this ecosystem. For 
example, members of the order Odonata are 
the major predator while the orders Diptera, 
Coleoptera and Lepidoptera are largely pests 
of the rice (Bambaradentya et al., 2004). In 
addition, aquatic insects could also have 
applications as _ bioindicators of the 
ecosystem (Federico Rizo-Patron et al., 
2013). Therefore, understanding taxonomy 
and biodiversity of these insects is important 
for applications as well as appropriate and 


INTRODUCTION 


Rice fields are an important artificial 
wetland ecosystem possessing different 
categories of the environmental conditions 
and providing habitats for different faunas 
(Lawler, 2001; Bambaradeniya and 
Amerasinghe, 2003). Diverse groups of 
organisms from microbial, fungi, plankton 
and macroinvertebrates to large vertebrates 
such as fishes, snakes and birds utilize this 
artificial ecosystem (Edirisinghe and 
Bambaradeniya, 2010). It has been found 
that biological diversity is positively related 
with the productivity of agricultural 
ecosystems (Roger et al., 1991; Scherr and 
McNeely, 2008; Tscharntke et al., 2012). 
However, using chemical pesticides and 





170 TROPICAL NATURAL HISTORY 20(2), AUGUST 2020 


sustainable management of the rice field 
ecosystem. 

The traditional method of species 
identification based on morphological 
characters has some limitations such as 
incomplete taxonomic knowledge, high 
intraspecific morphological variation, a lack 
of valid morphological characters for 
species recognition of particular life stages 
and damaged specimens. Beside these 
limitations, insufficient expertise or lack of 
an expert who has experience in the 
taxonomy of relevant insect groups also 
present major obstacles. One possible way 
to solve these issues is using the DNA 
barcode. Although there are also some 
limitations of using molecular genetic data 
for species identification (Frézal and 
Leblois, 2008) this method has more 
advantages | compared to traditional 
taxonomy (Sweeney et al., 2011, Gill et al., 
2013). DNA barcoding has been used 
successfully for species identification and 
documenting biodiversity of many insect 
groups (Rach et al., 2008; Damm et al., 
2010; Jinbo et al., 2011; Sweeney et al., 
2011; Jackson et al., 2014). Thus, this DNA 
marker can potentially be used for 
identification and documenting species 
diversity of rice field aquatic insects. 

Thailand is the major rice exporting 
country. Approximately 1896 (about 9 
million hectares) of the land is used for rice 
cultivation (http://www.ricethailand.go.th). 
However, data on biodiversity of the insects 
in rice fields is scanty due largely to the lack 
of the taxonomic knowledge. In this study, 
we used DNA barcoding sequences to test 
the efficiency of this molecular marker for 
species identification using public DNA 
barcode libraries (1.e. BOLD and GenBank). 
We also used the species delimitation 
approach based on the mitochondrial 
cytochrome c oxidase I (COI) sequence to 








document species diversity in the rice field 
ecosystem in northeastern Thailand. 


MATERIALS AND METHODS 


Specimen collection and identification 

Insects were collected from the eight rice 
fields in four provinces in northeastern 
Thailand between September and November 
2017 (Table 1). Specimens were collected 
using the kick-sampling method with a 500- 
um mesh hand net. Adult flying insects 
were collected using a sweep net along a 
line transect across the rice field. To 
compare the efficiency of DNA barcoding 
with morphological taxonomy in real 
practice (i.e. when monitoring / examining 
species diversity in the rice field, which 
usually involves people who are not an 
expert in particular insect groups), 
identifications of the specimens were carried 
out by university students who had skills at 
the taxonomic level of order / family of 
aquatic insects. Specimens were identified 
using the available keys and descriptions of 
the aquatic insects in Thailand (Sangpradub 
and Boonsoong, 2006; Boonsoong, 2014; 
Sangpradub, 2016). 
DNA extraction, polymerase chain reaction 
(PCR) and sequencing 

Genomic DNA was extracted from larva, 
nymph or soft tissue part of the adult using 
the Genomic DNA Extraction mini kit (GF- 
1 Tissue DNA Extraction Kit, Vivantis, 
Selangor Darul Ehsan, Malaysia). Polymerase 
chain reactions (PCR) were conducted for 
cytochrome c oxidase subunit I (COD using 
primers LCO1490 (5'-GGTCAACAAATC 
ATAAAGATATTGG-3’) and HCO2198 (5'- 
TAAACTTCAGGGTGACCAAAAAATCA- 
3^) (Folmer et al., 1994). PCR reaction and 
conditions followed the method of Rivera 
and Currie (2009). Agarose gel electrophoresis 
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TABLE 1. Sampling locations of the aquatic insects used in this study 


Location —  — — Code  Latitude/Longitude —— Date ` 
Yang Ta Lad, Kalasin Province KSI 16? 29' 13" N/103? 12' 42" E 30 Sep 2017 
Khao Wong 1, Kalasin Province KS2 16? 38' 22" N/104? 08' 39" E 07 Oct 2017 
Khao Wong 2, Kalasin Province KS3 16° 32' 09" N/104? 07' 15" E 07 Oct 2017 
Yang Sisurat 1, Maha Sarakham Province MKI 15? 40' 23" N/103? 09' 48" E 28 Oct 2017 
Yang Sisurat 2, Maha Sarakham Province MK2 15? 36' 15" N/103? 09' 41" E 28 Oct 2017 
Na Dun, Maha Sarakham Province MK3 15? 41' 00" N/103? 14' 10" E 28 Oct 2017 
Mueang Amnat Charoen, Amnat Charoen Province AM 15? 49' 03" N/104? 37' 48" E 4 Nov 2017 
Prangku, Si Sa Ket Province SK 14? 48' 59" N/104? 04' 00" E 23 Oct 2017 


was used to check the PCR products. The 
HiYield Gel/PCR DNA Extraction Kit 
(RBC BIOSCIENCE, Taiwan) was used to 
purified the PCR products and send for 
sequencing at 1‘ Base sequencing service 
(Malaysia) using the same primers as in the 
PCR. 
Data analysis 

A total of 98 sequences were included in 
data analyses. Sequences were deposited in 
GenBank under accession numbers: 
MH$881228-MH881247 and MHS$881249- 
MH881326. To assess the species identification 
based on COI DNA sequence data, the 
"identification" menu with *Species Level 
Barcode Records" option in Barcoding of 
Life Data Systems (http://v4.boldsystems.org) 
was used. In addition, we also used the 
Poisson tree processes (PTP) (Zhang et al., 
2013) to delimit species of rice field insects 
and to investigate hidden diversity in 
morphological species. The maximum 
likelihood tree was inferred in the RAxML 
web server version (https://embnet.vital- 
it.ch/raxml-bb/) (Stamatakis et al., 2008). 
The output tree was used in PTP species 
delimitation analyses using the web server 
version (https://species.h-its.org/) (Zhang et 
al., 2013). The PTP analysis was run for 
100,000 MCMC generations using a 

















thinning value of 100 and 10% burn-in. In 
addition to the tree-based (1.e. PTP) we also 
used a distance-based species delimitation 
method called Automatic Barcode Gap 
Discovery (ABGD) (Puillandre et al., 2012). 
ABGD analysis was performed in web 
server version (https://bioinfo.mnhn.fr/abi/ 
public/abgd/abgdweb.html). Criteria setting 
for ABGD analysis as follow: prior range 
for maximum intraspecific divergence 
(0.001, 0.1), minimum slope increase (X) of 
1.0 and genetic distance calculated based on 
K80 corrected distances. 

To indicate phylogenetic relationships 
between sequences of the rice field insects 
included in this study, three methods of 
phylogenetic analyses were conducted 
including the neighbor-joining (NJ) tree, the 
maximum likelihood (ML) and the Bayesian 
analysis (BA). Because specimens were 
from distantly related taxa, the phylogenetic 
trees were analyzed separately for 
specimens from each order (i.e. Odonata, 
Ephemeroptera and Hemiptera). The NJ 
analysis was performed in the MEGA 6 
(Tamura et al., 2013) based on Kimura 2- 
parameter (K2P) model. The branch supports 
were calculated based on bootstrapping 
method with 1,000 replications. The ML 
tree was estimated in PhyML 3.0 (Guindon 
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et al., 2010). Branch support was calculated 
using approximate likelihood ratio tests 
(Anisimova and Gascuel, 2006). Bayesian 
analysis was performed with MrBayes 3.04b 
(Huelsenbeck and Ronquist, 2001), and was 
run for 2,000,000 generations with sampling 
every of 100 generations and the first 25% 
were discarded as burn-in. For ML and BA 
analyses, the jModelTest (Darriba et al., 
2012) was used to select the best-fit 
substitution model for sequences data based 
on the corrected Akaike information 
criterion. The COI sequences of Callibaetis 
ferrugineus (GenBank accession number 
KT707077), Procladius sp. and Tipula sp. 
were used as outgroup for phylogenetic 
analyses of Odonata, Hemiptera and 
Ephemeroptera, respectively. 


RESULTS 


Morphological identification 

A total of 175 specimens were 
morphologically examined and 98 were 
used for molecular analysis (Table 2). Based 
on morphological characters, only 32 
specimens were successfully identified to 
species level (Table 2). All morphological 
identifications that could assign a specimen 
to a species were for the adults of the order 
Odonata, excepting one adult specimen of 
Hemiptera that was identified as Diplonychus 
rusticus Fabricius (Table 2). In total, 13 
species were recognized. Seven specimens 
were assigned into genus and the remainder 
could only be assigned to order or family 
levels. 
DNA barcode species identification and 
species delimitation 

A total of 98 specimens were included in 
DNA barcode analysis (Table 2). 
Identification of the specimens based on 
DNA barcode revealed that 63 of 98 














sequences were successfully identified into 
species level (Table 2). These sequences 
were belonging to 13 insect species, most 
(62 sequences) of these were from the order 
Odonata. For other insect orders (12 
specimens of Hemiptera, 3 specimens of 
Diptera and 5 specimens of Ephemeroptera), 
only single sequence of the order Diptera 
was successfully identified in the BOLD 
(Table 2). The results of DNA barcode and 
species delimitation of each insect order 
based on PTP and ABGD methods are as 
follows.— 
Odonata 

A total of 78 sequences were obtained 
and were morphologically identified into 11 
species. DNA barcoding in BOLD revealed 
that 62 specimens were successfully 
identified into 12 species and 16 specimens 
were ambiguous identifications (Table 2). 
The NJ tree (Fig. 1) revealed that all species 
were formed monophyletic group agree with 
DNA barcode species identification. The 
PTP and ABGD species delimitation 
method revealed that the samples included 
17 and 16 species, respectively (Fig. 1). Two 
specimens of Pseudagrion australasiae 
Sélys were delimited into two species 
according to the PTP analysis (Fig 1) but 
ABGD combined them into single species. 
Ephemeroptera 

Five sequences from the larvae of 
Ephemeroptera were obtained. These 
specimens are only morphologically identified 
to the level of order. DNA  barcoding 
identification of these sequences in BOLD 
were not successful because no match 
sequences were found. The NJ tree (Fig. 2) 
revealed that all five specimens were 
separated into two groups. Four specimens 
formed one group and remaining formed 
another. Both PTP and ABGD species 
delimitations revealed that these groups 
represent different species (Fig. 2). 
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TABLE 2. Specimens used for DNA barcode analysis and the species identification based on morphology and 








DNA barcoding sequence 
Location Specimen GenBank 
code code accession 
number 

AM 0do1316 MH881259 
AM 0do1317 MH881260 
AM odo1312 MH881271 
AM 0do1319 MH881318 
AM odol311 MH881288 
AM 0do132 MH881289 
AM 0do1321 MH881287 
AM 0do1325 MH881290 
AM 0do1326 MH881291 
AM 0do1331 MH881292 
AM 0do131 MH881293 
AM 0do1329 MH881294 
AM 0do1327 MH881295 
AM 0do135 MH881296 
AM 0do137 MH881297 
AM 0do138 MH881286 
AM 0do139 MH881298 
AM 0do1314 MH881326 
AM 0do134 MH881249 
AM 0do136 MH881250 
AM 0do1332 MH881252 
AM 0do1313 MH881254 
AM 0do1330 MH881253 
AM 0do1320 MH881251 
KS1 2APAI MH881256 
KSI 0do037 MH881270 
KSI Cul01 MH881247 
KSI Hemi01 MH881242 
KSI Hele04 MH881243 
KSI Noto012 MH881234 
KSI Noto022 MH881235 
KSI Noto01 MH881236 
KSI Noto013 MH881237 
KS1 Tipu03 MH881245 
KS2 0do53 MH881302 
KS2 odocal62 MH881303 
KS2 odo55 MH881281 
KS2 odo52 MH881282 
KS2 1CIN1 MH881304 


Life stage 


Larva 
Larva 
Larva 
Larva 


Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 


Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Adult 


Larva 
Larva 
Adult 
Adult 
Nymph 
Nymph 
Nymph 


Nymph 


Larva 
Larva 
Larva 
Larva 


Larva 


Adult 


Morphological 
identification 


Odonata 
Odonata 
Odonata 
Odonata 


Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 
Odonata 


Odonata 

Odonata 

Odonata 

Odonata 

Odonata 

Odonata 

Odonata, 
Libellulidae, Acisoma 
panorpoides 
Odonata 

Diptera, Culicidae, 
Culex sp. 

Hemiptera 
Hemiptera, 
Naucoridae, 
Heleocoris sp. 
Hemiptera, 
Notonectidae, Nychia 
sp. 

Hemiptera, 
Notonectidae, Nychia 
sp. 

Hemiptera, 
Notonectidae, Nychia 
sp. 

Hemiptera, 
Notonectidae, Nychia 
sp. 

Diptera, Tipulidae 
Odonata 

Odonata 

Odonata 


Odonata 


Odonata, 
Coenagrionidae 
Ceriagrion 
cerinorubellum 


DNA barcode 
identification 


Acisoma panorpoides 
Acisoma panorpoides 
Agriocnemis minima 
Crocothemis servilia 
Pantala sp. 

Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Ischnura senegalensis 
Pseudagrion 
microcephalum 
Pseudagrion rubriceps 
Pseudagrion rubriceps 
Pseudagrion rubriceps 
Pseudagrion rubriceps 
Pseudagrion rubriceps 
Pseudagrion rubriceps 
Acisoma panorpoides 


Agriocnemis minima 
Culex bitaeniorhynchus 


No match 
No match 


No match 


No match 


No match 


No match 


No match 

Aciagrion pallidum 
Aciagrion pallidum 
Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 
Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 
Ceriagrion 
cerinorubellum 


Percentage 
closest 
match in 
BOLD 
99.29 

99.47 

99.46 

99.13 


100 
100 
98.30 
100 
100 
100 
100 
100 
99.82 
99.56 
100 
99.15 
99.56 
99.82 


99.63 
99.63 
98.56 
99.46 
100 

99.28 
99.11 


98.92 
99.11 


174 
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TABLE 2. Continue. 


Location 
code 


KS2 


KS2 


KS2 


KS2 


KS2 


KS2 


KS2 


KS2 


KS2 
KS2 
KS2 
KS2 
KS3 


KS3 


KS3 


KS3 


KS3 


KS3 


KS3 


KS3 


KS3 


KS3 


KS3 


MKI 


MKI 


Specimen GenBank Life stage Morphological DNA barcode 
code accession identification identification 
number 
ICSEMI MH881316 Adult Odonata, Libellulidae Crocothemis servilia 
Crocothemis servilia Pantala sp. 
1CSE2FM MH881317 Adult Odonata, Libellulidae Crocothemis servilia 
Crocothemis servilia Pantala sp. 
1DNE3 MH881306 Adult Odonata, Diplacodes nebulosa 
Libellulidae, 
Diplacodes nebulosa 
1DNE2 MH881305 Adult Odonata, Diplacodes nebulosa 
Libellulidae, 
Diplacodes nebulosa 
1ISE1 MH881299 Adult Odonata, Ischnura senegalensis 
Coenagrionidae, 
Ischnura 
senegalensis 
1ISE2 MH881300 Adult Odonata, Ischnura senegalensis 
Coenagrionidae, 
Ischnura 
senegalensis 
INTUI MH$881320 Adult Odonata, Neurothemis intermedia 
Libellulidae, 
Neurothemis tullia 
INTU2 MH881321 Adult Odonata, Neurothemis intermedia 
Libellulidae, 
Neurothemis tullia 
Eph063 MH881229 Larva Ephemeroptera No match 
Eph62 MH881230 Larva Ephemeroptera No match 
Eph51 MH881231 Larva Ephemeroptera No match 
Eph52 MH881232 Larva Ephemeroptera No match 
1APY3 MH881272 Adult Odonata, Agriocnemis minima 
Coenagrionidae, 
Agriocnemis minima 
odo78 MH881279 Larva Odonata Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 
odo79 MH881280 Larva Odonata Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 
odo82 MH881283 Larva Odonata Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 
ICOLI MH881315 Adult Odonata, Ceriagrion 
Coenagrionidae, coromandelianum 
Ceriagrion 
coromandelianum 
INTUM2 MH881325 Adult Odonata, Neurothemis tullia 
Libellulidae, 
Neurothemis tullia 
Rus08 MH881244 Adult Hemiptera, No match 
Belostomatidae, 
Diplonychus rusticus 
Tipo08 MH881246 Larva Diptera, No match 
Chironomidae 
Noto) MH881233 Nymph Hemiptera, No match 
Notonectidae, 
Nychia sp. 
odo83 MH881284 Larva Odonata Pseudagrion australasiae 
Pseudagrion 
malabaricum 
odo84 MH881285 Larva Odonata Pseudagrion australasiae 
Pseudagrion 
malabaricum 
2APA2 MH881257 Adult Odonata, Acisoma panorpoides 
Libellulidae, 
Acisoma panorpoides 
odo91 MH881255 Larva Odonata Acisoma panorpoides 


Percentage 
closest 
match in 
BOLD 
99.65 


99.47 


98.43 


98.84 


99.82 


99.56 


100 


98.67 


100 


99,25 


98.40 


99.43 


98.57 


98.93 
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TABLE 2. Continue. 


Location 
code 


MKI 
MKI 
MKI 
MKI 
MKI 
MKI 
MKI 
MKI 
MKI 
MKI 





MKI 


MKI 
MK2 


MK2 


MK2 


MK2 


MK2 


MK2 


MK2 


MK2 


MK3 


MK3 


MK3 
MK3 
MK3 
MK3 
MK3 
SK 


SK 


SK 


Specimen 
code 


odo99 
odo910 
odo92 
odo93 
odo94 
odo95 
odo96 
0do97 
odo98 
3APYII 


2DTR2 


Eph91 
3APYI 


3DNE3 


3DTRI 


3DTR2 


3DTR3 


3DTR4 


2DTRI 


GPR 


odol 13 


odol12 


odol 14 
Hemil11 
Hemi112 
Hemil13 
Hemi1 14 
4APAI 


4APYI 


4APY2 


GenBank 
accession 
number 


MH881261 
MH881262 
MH881263 
MH88 1264 
MH881265 
MH88 1266 
MH88 1267 
MH881268 
MH88 1269 
MH881273 





MH881312 


MH881228 
MH881276 


MH88 1307 


MH88 1308 


MH88 1309 


MH881310 


MH881311 


MH881314 


MH881319 


MH881277 


MH881278 


MH881301 
MH881238 
MH881239 
MH881240 
MH881241 
MH881258 





MH88 1274 


MH881275 


Life stage 


Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Larva 
Adult 


Adult 


Larva 
Adult 


Adult 


Adult 


Adult 


Adult 


Adult 


Adult 


Adult 


Larva 


Larva 


Larva 
Adult 
Adult 
Adult 
Adult 
Adult 


Adult 


Adult 


Morphological 
identification 


Odonata 

Odonata 

Odonata 

Odonata 

Odonata 

Odonata 

Odonata 

Odonata 

Odonata 

Odonata, 
Coenagrionidae, 
Agriocnemis 
pygmaea 

Odonata, 
Libellulidae, 
Diplacodes trivialis 
Ephemeroptera 
Odonata, 
Coenagrionidae, 
Agriocnemis 
pygmaea 

Odonata, 
Libellulidae, 
Diplacodes nebulosa 
Odonata, 
Libellulidae, 
Diplacodes trivialis 
Odonata, 
Libellulidae, 
Diplacodes trivialis 
Odonata, 
Libellulidae, 
Diplacodes trivialis 
Odonata, 
Libellulidae, 
Diplacodes trivialis 
Odonata, 
Libellulidae, 
Diplacodes trivialis 
Odonata, 
Libellulidae, 
Orthetrum pruinosum 


Odonata 


Odonata 


Odonata 
Hemiptera 
Hemiptera 
Hemiptera 
Hemiptera 
Odonata, 
Libellulidae, Acisoma 
panorpoides 
Odonata, 
Coenagrionidae, 
Agriocnemis 
pygmaea 
Odonata, 
Coenagrionidae, 
Agriocnemis 
pygmaea 


DNA barcode 
identification 


Agriocnemis minima 
Agriocnemis minima 
Agriocnemis minima 
Agriocnemis minima 
Agriocnemis minima 
Agriocnemis minima 
Agriocnemis minima 
Agriocnemis minima 
Agriocnemis minima 
Agriocnemis pygmaea 


Diplacodes trivialis 


No match 
Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 


Diplacodes nebulosa 


Diplacodes trivialis 


Diplacodes trivialis 


Diplacodes trivialis 


Diplacodes trivialis 


Diplacodes trivialis 


Orthetrum pruinosum 
Orthetrum pruinosum 
neglectum 

Orthetrum testaceum 
Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 
Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 
Ischnura senegalensis 
No match 

No match 

No match 

No match 

Acisoma panorpoides 


Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 


Agriocnemis pygmaea 
Ischnura senegalensis 
Agriocnemis femina 


Percentage 
closest 
match in 
BOLD 


99.28 
98.91 
99.09 
98.01 
99.82 
99.64 
99.64 
99.82 
99.46 
100 


100 


99.22 


99.61 


99.03 


99.40 


99.2] 


99.79 


99.81 


176 
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Location Specimen GenBank Life stage Morphological DNA barcode Percentage 
code code accession identification identification closest 
number match in 
BOLD 
SK 4DTRI MH881313 Adult Odonata, Diplacodes trivialis 99.40 
Libellulidae, 
Diplacodes trivialis 
SK 4NTU21 MH881322 Adult Odonata, Neurothemis tullia 98.57 
Libellulidae, 
Neurothemis tullia 
SK 4NTU3 MH88 1323 Adult Odonata, Neurothemis tullia 98.04 
Libellulidae, 
Neurothemis tullia 
SK 4NTU31 MH88 1324 Adult Odonata, Neurothemis tullia 98.22 
Libellulidae, 


Neurothemis tullia 


Details of location code are provided in Table 1. 





Hemiptera 

Twelve sequences were obtained from 
the order Hemiptera. Morphological 
identifications revealed that they are 


comprised of D. rusticus (n-l), genus 
Nychia Stal (n=5) and Heleocoris Stal (n=1) 
and remaining five specimens were only 
identified at the order level (Table 2). Like 
those of the Ephemeroptera, DNA barcode 
identification in BOLD found no match for 
sequences obtained from this order. NJ tree 
(Fig. 3) revealed three groups, all with 
strong supported. These groups comprised 
of four species based on PTP species 
delimitation analysis. However, the ABGD 
method delimits specimens of the order 
Hemiptera into three species (Fig. 3). All 
five sequences from the specimens that were 
morphologically identified at the order level 
formed one species. The genus Nychia and 
genus Heleocortis formed two separated 
species with D. rusticus forming another 
(Fig. 3). 
Diptera 

Only three sequences were obtained 
from the order Diptera. Two specimens 
were morphologically identified into genus 
Tipula Linnaeus of the family Tipulidae and 
one specimen was identified into the genus 

















Culex Linnaeus of the family Culicidae 
(Table 2). DNA barcode analysis in BOLD 
found that the Culex specimen was match 
with Culex bitaeniorhychus Giles with high 
percent similarity (99.11%). However, two 
sequences of the genus Tipula found no 
match in BOLD. The PTP species 
delimitation revealed two but the ABGD 
method found one species of Diptera. Two 
sequences of the genus Tipula belonged to 
the same species and the Culex specimen 
formed another species. 


DISCUSSION 


Aquatic insects play diverse roles in the 
rice field ecosystem. Therefore, knowledge 
of taxonomic and species diversity is crucial 
for all aspects of further study of these 
organisms. We have revealed that 
mitochondrial COI sequences could be 
successfully used for species identification 
of one of the major groups of rice field 
insects, the order Odonata. The major 
advantage of DNA barcoding is that it can 
be used to associate different life stages 
(Miller et al., 2005; Gattolliat and Monaghan, 
2010; Pramual and Wongpakam, 2014). 
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odo93 Agriocnemis minima 




















odo97 Agriocnemis minima 

odo95 Agriocnemis minima 

odo99 Agriocnemis minima 

o0do910 Agriocnemis minima 

0do037 Agriocnemis minima Agriocnemis minima 
odo96 Agriocnemis minima 
ee odo98 Agriocnemis minima 
0do1312 Agriocnemis minima 
1APY3 Agriocnemis minima 
0do112 Agriocnemis pygmaea 
100/100100 3APY1 Agriocnemis pygmaea 

odo113 Agriocnemis pygmaea 

0do78 Agriocnemis pygmaea Agriocnemis pygmaea 
4APY2 Agriocnemis pygmaea 

4APY1 Agriocnemis pygmaea 

odo79 Agriocnemis pygmaea 


0do1314 Pseudagrion microcephalum I Pseudagrion microcephalum 





odo83 Pseudagrion australasiae Pseudagrion australasiae # 1 


odo84 Pseudagrion australasiae Pseudagrion australasiae # 2 


P. australasiae 


0do1332 Pseudagrion rubriceps 
100/100/1.00 
93/94/0.99 


0do136 Pseudagrion rubriceps 

0do1330 Pseudagrion rubriceps 

0do1320 Pseudagrion rubriceps NiseHdggronmmcops 
0do134 Pseudagrion rubriceps 


0do1313 Pseudagrion rubriceps 


VE 1CIN1 Ceriagrion cerinorubellum BE Ceriagrion cerinorubellum 


1COL1 Ceriagrion coromandelianum 8 W Ceriagrion coromandelianum 
odo53 Aciagrion pallidum 
100/100/1:00 Vm . | Aciagrion pallidum 
odocal62 Aciagrion pallidum 
89/100/1.00 r- 0d0138 /schnura senegalensis 
0do1321 Ischnura senegalensis 


100/100/1:00 0do1311 Ischnura senegalensis 


Ischnura senegalensis 
0do132 Ischnura senegalensis 
0do135 Ischnura senegalensis 
1ISE1 Ischnura senegalensis 
0do1319 Crocothemis servilia 
100/99/1.00 
1CSEM1 Crocothemis servilia Crocothemis servilia 
1CSE2FM Crocothemis servilia 
77/98/0.974NTU3 Neurothemis tullia 
4NTU31 Neurothemis tullia 
100/100/1.00 L 4NTU21 Neurothemis tullia 


100/98/1.00 1NTUM?2 Neurothemis tullia 


Neurothemis tullia 


100/95/1.00 1NTU1 Neurothemis intermedia SG 
Neurothemis intermedia 
1NTU2 Neurothemis intermedia 
GPR Orthetrum pruinosum GR Orthetrum pruinosum 
1DNE2 Diplacodes nebulosa 
100/100/1.00 |- 3DNE3 Diplacodes nebulosa Diplacodes nebulosa 
1DNE3 Diplacodes nebulosa 
3DTR1 Diplacodes trivialis 
100/100/1.00 | 3 TR? Diplacodes trivialis 
93/93/0.73| r 3DTR3 Diplacodes trivialis . ux 
. m Diplacodes trivialis 
2DTR2 Diplacodes trivialis 
98/90/1.00] | 4DTR1 Diplacodes trivialis 
97/91/0.83_. 2DTR1 Diplacodes trivialis 


2APA2 Acisoma panorpoides 





100/87/93 





2APA1 Acisoma panorpoides 


4APA1 Acisoma panorpoides 


N p Acisoma panorpoides 
odo91 Acisoma panorpoides 


0do1316 Acisoma panorpoides 
81/98/0.96 


0do1317 Acisoma panorpoides 





KT707077 Callibaetis ferrugineus ll Outgroup 


a 
0.05 


FIGURE 1. Neighbor-joining tree for 61 haplotypes of the mitochondrial cytochrome c oxidase I (COD 
sequences of Odonata collected from rice fields, Thailand. Bootstrap values for neighbor-joining and maximum 
likelihood (ML) and posterior probability of Bayesian analysis are shown above or near the branches. Scale bar 
represents 0.05 substitutions per nucleotide position. Vertical bar indicates species delimited by Poisson tree 
processes (PTP) (grey) and Automatic Barcode Gap Discovery (ABGD) (black) methods. 
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81/81/0.94| Eph62 Ephemeroptera 
88/85/0.97 








Eph52 Ephemeroptera 

Eph51 Ephemeroptera Species 1 

Eph063 Ephemeroptera 

Eph91 Ephemeroptera "RH Species 2 


Tipu03 Tipula sp. BOutgroup 


0.02 


FIGURE 2. Neighbor-joining tree for 5 sequences of the mitochondrial cytochrome c oxidase I (COI) 
sequences of 4 taxa of Ephemeroptera from Thailand. Bootstrap values for neighbor-joining and maximum 
likelihood (ML) and posterior probability of Bayesian analysis are shown above or near the branches. Scale 
bar represents 0.05 substitutions per nucleotide position. Vertical bar indicates species delimited by Poisson 
tree processes (PTP) (grey) and Automatic Barcode Gap Discovery (ABGD) (black) methods. 


Hemi111 Hemiptera 






Hemi114 Hemiptera 


Hemi113 Hemiptera Species 1 


100/100/1.00 






Hemi112 Hemiptera 


95/96/1.00 Hemi01 Hemiptera 










100/100/1.00 Hele04 Heleocoris sp. E Heleocoris sp. 







Rus08 Diplonychus rusticus [ species 2 


Diplonychus rusticus 


Noto81 Nychia sp. 










100/95/1.00 
Noto01 Nychia sp. 


Noto012 Nychia sp. 
Noto022 Nychia sp. 
Noto013 Nychia sp. 


Nychia sp. 


TipoO8 Procladius sp. B Outgroup 


0.02 


FIGURE 3. Neighbor-joining tree for 12 sequences of the mitochondrial cytochrome c oxidase I (COI) 
sequences of two taxa of Hemiptera from Thailand. Bootstrap values for neighbor-joining and maximum 
likelihood (ML) and posterior probability of Bayesian analysis are shown above or near the branches. Scale 
bar represents 0.05 substitutions per nucleotide position. Vertical bar indicates species delimited by Poisson 
tree processes (PTP) (grey) and Automatic Barcode Gap Discovery (ABGD) (black) methods. 


This is highly useful for the insects which larvae with three conspecific adults. 


are taxonomically challenging in some life 
stages. For example, identification of the 
larva of Odonata based on morphology is 
very difficult and all of our specimens could 
only be assigned to the level of order. In 
contrast, DNA barcodes can associate 33 


Although the remaining (9) larvae were not 
associated with our adult sequences, they 
could be successfully identified to species in 
DNA barcode database (1.e. BOLD). 

Despite a considerable rate (8096, 62 
sequences) of successful identification for 
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Odonata in BOLD, some (20%, 16 
sequences) were ambiguous and raised the 
problem of the uncertainty of the taxonomic 
identity of reference sequences in the public 
database (Pramual et al, 2016). These 
sequences are matching with more than one 
species in the database. It could be possible 
that some of these species (Agriocnemis 
femina Brauer and A. pygmaea Rambur; P. 
australasiae and P. malabaricum Fraser) are 
closely related species (O'Grady and May, 
2003) thus, lineage sorting might not yet be 
complete (Funk and  Omland, 2003) 
resulting in ambiguous identification. of 
these specimens. Closely related species are 
often problematic for DNA barcoding of 
insects (Elias et al, 2007, Pramual and 
Adler, 2014) including members of the 
Odonata (Rach et al, 2008). Another 
possibility is that the sources of some 
sequences deposited in the public database 
(e.g. BOLD and GenBank) might be 
incorrectly identified. For example, our 
specimens of Crocothemis servilia Drury 
were match with C. servilia and Pantala sp. 
in BOLD. These two genera could clearly 
be differentiated based on morphological 
characters. Thus, ambiguous identification 
of these specimens is most likely due to 
erroneous identification of the sources of 
sequences deposited in the public database. 
This issue has already been raised by several 
studies (e.g. Nilsson et al., 2006, Shen et al., 
2013; Pramual et al., 2016) and suggests 
some limitations of DNA barcoding. 
Specimens from other orders including 
Diptera, Ephemeroptera and  Hemiptera 
were far less successfully identified in 
BOLD; just one sequence from Diptera (C. 
bitaeniorhychus) was successfully identified. 
This is because of a lack of reference 
sequences in the database. One possible way 
of documenting biodiversity if there are no 
reference sequences is to use the DNA- 


























based, species delimitation method. The 
PTP and ABGD analyses revealed 25 and 
23 species, respectively, among 98 
sequences included in the analysis. There is 
one more species from the Coleoptera 
included in this study but PTP analysis 
could not be performed because only a 
single specimen was included in this study. 
Thus, in total we found 26 of the aquatic 
insect species in the rice field among 99 
specimens. Although this number of species 
is far less than the possible diversity of 
insects in the rice field ecosystem, it was 
more than twice that based on 
morphological identification (13 species). 
This indicates that DNA barcoding 
sequences are potentially useful for 
documenting biodiversity in this under 
investigated ecosystem in Thailand. 

In conclusion, we found that DNA 
barcode sequences are successful for 
identification of Odonata in the rice field as 
several species have reference sequences in 
the public database (e.g. BOLD). However, 
the approach had far less success with other 
insect orders because of limitations of the 
reference sequence library. Therefore, this 
should be a priority topic for study in the 
future. When the species are unknown, we 
have demonstrated in this study that 
numbers of species can still be documented 
using the molecular species delimiting 
approach. Thus, despite lacking names for 
species, DNA barcode can be potentially 
used for biodiversity assessment. 
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